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Abstract 

The N-body problem for one-dimensional self-gravitating systems 
has been often studied to test theories of galactic evolution and sta- 
tistical mechanics. We consider the general relativistic version of this 
system and obtain the first exact solution to the 2-body problem in 
which spacetime is not flat. In the equal mass case we obtain an 
explicit expression for the proper separation of the two masses as a 
function of their mutual proper time. 
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iV-body self-gravitating systems have a long history in physics and are of 
interest in studying both star systems (small N > 2) and galactic evolution 
(large N). One- dimensional (or lineal) models of such systems have been 
of particular interest |]] in that they avoid some difficulties due to three 
dimensions, including evaporation, singularities, and energy dissipation in the 
form of gravitational radiation, as well as admitting a level of computational 
and analytic analysis which is dramatically simpler. 

We consider in this paper the iV-body problem for a relativistic self- 
gravitating lineal system and formulate a canonical theory determining its 
Hamiltonian. In the 2-body case we obtain an exact solution (valid to all 
orders in the gravitational coupling k) in which spacetime outside the mov- 
ing matter sources is not flat. In the equal mass case we obtain an explicit 
expression for the proper separation of the two point masses as a function 
of their mutual proper time. This is the first non-perturbative relativistic 
curved-spacetime treatment of this problem, providing new avenues for in- 
vestigation of one- dimensional self-gravitating systems. 

For our lineal self-gravitating system we choose a modification of Jackiw- 
Teitelboim lineal gravity, in which the scalar curvature is equated to a cos- 
mological constant 

R-A = (1) 

in the absence of other matter fields |J . This model has been of considerable 
interest as a model theory of quantum gravity 0. Here we couple iV point 
masses to this theory so that we have a generally covariant self-gravitating 
system with non-zero curvature outside the point sources. We do not include 
collisional terms, so that the bodies pass through each other. 

Since the Einstein action is a topological invariant in 2 spacetime dimen- 
sions, we must incorporate a scalar (dilaton) field into our action, which we 
take to be 



/ d2x {* r + ^v^v^ + a} 



+ Ip (2) 



where ^ is the dilaton field, and g are the metric and its determinant and 
R is the Ricci scalar, with k = 87rG/c 4 . Ip is the action of A" point masses 
minimally coupled to gravity 

Ip = - J d 2 x ^m a J dr a l-g^{x) d ** d ** \ 5 2 {x - z a (r a )) 



a=l 



1 



where r a is the proper time of a-th particle. Variation of the action (0) with 
respect to the metric, dilaton field, and particle coordinates yields the field 
equations 



~V M *V^ - g, v Q V A v&V A vI/ - V 2 *) - V M V,* = k7£ + | 
where the stress-energy due to the point masses is 



N f 1 cfc°" <fc p 

T £ = E "W dr a -=g^g up -^-^6 2 (x - z a (r a )) (5) 



and is conserved. We observe that (|3]) is a closed system of N + 1 equations 
for which one can solve for the single metric degree of freedom and the N 
degrees of freedom of the point masses; it reduces to ([I]) if all masses vanish. 
The evolution of the dilaton field is governed by the evolution of the point- 
masses via (|]). The left-hand side of (f|) is divergenceless (consistent with the 
conservation of T^), yielding only one independent equation to determine 
the single degree of freedom of the dilaton. 

Working in the canonical formalism we make use of the decomposition 
y/^R = — 2<9 ( v /7-K') + 2di(^N 1 K — 7 -1 <9iiV ) where the extrinsic curva- 
ture K = (2N 'y)~ 1 (2diN 1 — 7~ 1 iV 1 <9i7 — <9 7), an d rewrite the action (0) in 
the form |4]] 

I = Jdx 2 !^p a zj(x - z a (x )) + 7T 7 + + N R° + iVii? 1 j (6) 

where 7 = gn,N = (— g 00 ) -1 / 2 , iVi = g w , n and II are conjugate momenta 
to 7 and \I> respectively. The quantities No and N\ are Lagrange multipliers 
which enforce the constraints R° = = R l , where 



(\b>) 2 ( \fj> \' A [75 
R° = -Ky/^ + 2Ky/^U + j-^-- — ) + ^-J2J^ + ml5(x-z a (x )) 

R 1 = ^7r--n*' + 2tt' + y^5(x-z a (x )) (7) 
7 7 a 7 



2 



with the symbols ( ' ) and ( ' ) denoting do and d\, respectively. We identify 
the dynamical and gauge (i.e. co-ordinate) degrees of freedom by writing 
the generator arising from the variation of the action at the boundaries in 
terms of (M/'/^/t")' and n', which we can easily solve for since these are the 
only linear terms in the constraints. We then fix the frame of the physical 
space-time coordinates in a manner similar to the (3 + l)-dimensional case 

1- 

Carrying out this procedure, we find that we can consistently choose the 
coordinate conditions 7 = 1 and II = 0. Eliminating the constraints, the 
action (|6]) then reduces to 

I = J d 2 x ^J2p a z a S(x - z a ) - 7i| . (8) 

where the reduced Hamiltonian is H = J dxH = — - J dxA^fr , where A = 
d 2 /dx 2 , and \I> = ^(x, z a , p a ) and is understood to be determined from the 
constraint equations which are now 

AM/ - VOT + «V -± + K J2<Jtt + ml5(x - z a ) = 

2V + 5> <f(z - *«) = 0. (9) 

a 

When A = the Hamiltonian reduces to that considered in ref. [[y] in the 
non-relativistic limit Q. 

For iV = 2 we solve these equations exactly following the method given 
in ref. 0. First, for z 2 < z%, we divide spacetime into three regions: Z\ < x 
((+) region), z 2 < x < Z\ ((0) region) and x < Z2 ((-) region) and set 
\& = — 41og|0| and 7r = x'- I n each region \ is a sum of terms linear in 
x and cf) obeys a harmonic oscillator equation, whose solutions are linear 
combinations of exponential functions of x in the ± regions and of either 
trigonmetric or exponential functions of x in the region depending on the 
size of A relative to the other integration constants. Matching these solutions 
at the boundaries x = z\ and x = z 2 of each region allows a determination of 
the coefficients of these linear combinations in the + and — regions in terms 
of those in the region. 

The magnitudes of both <ft and \ increase with increasing \x\ and so 
we must employ a boundary condition which guarantees that the surface 
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terms obtained in passing from (0) to vanish and the Hamiltonian remain 
finite. A perturbative analysis indicates that this condition must be \I/ 2 — 
4k, 2 x 2 + 2Ax 2 = C±x in the ± regions. Incorporating this into the matching 
conditions determines the coefficients of the linear combination in the region 
in terms of the momenta and positions of the bodies. 

This procedure is sufficient to solve for N , Ni and all other field variables 
exactly. Repeating the analysis for z\ < z 2 yields a similar solution with pi — ► 
— pi, and so the full solution is a function of (zi,pi), where pi = piSgn(zi — z?). 
The resultant expressions are rather long and cumbersome. However it is 
straightforward to show that the canonical equations imply conservation of 
the total momenta p\ + p 2 - Choosing a center of inertia frame with p\ = 
—P2 = p, yields considerable simplification of the exact solution, which we 



write employing the notation K = J (kX) 2 — ^, K — JYq — ^ and 









We obtain 



(f) ± (x) = A±e 



** Kx 4> (x) = A e^ Ko + B e~% 



1^0 



(10) 



in the (+), ( 



) and (0) regions respectively, where 




(11) 



4^o 



1/2 



{^0,^0} 



e 4 



The metric components are 



N 0{±) 



A(j> 2 ± N 0{0) = A<p 2 {0) 




(12) 



A2 K x _ tYI e -K x 



where A = 16 ]*° x e ^ zl ~ Z2 \ D = nf§ and 

'± = 2 (I + f ) <*' ± (I - f ) Gar ± ^)-|^^(— ' 

(13) 

The solution for z\ < z 2 is obtained by interchanging the suffices 1 and 2 
in the preceding solution. 

The matching conditions at (zi,z 2 ) force the relation 



K X K 2 = MiM 2 e K ^ Zl - Z2 \ (14) 
which determines X. From this the Hamiltonian 

H = -±(dx^ = -Uvr oo = — ■ (15) 



may be explicitly determined as a function of the coordinates and momenta 
of the particles. 

The parameter e = ±1, and is a constant of integration associated with 
the metric degree of freedom. Under time reversal, solutions with e = 1 
transform into those with e = — 1, ensuring invariance of the whole theory 
under this symmetry. It is straightforward to show from this solution that 
the Ricci scalar is equal to the cosmological constant everywhere except at 
the locations (zi(t), z 2 (t)) of the point masses. Hamilton's equations imply 
that 

AX K Q Ky K 2 

P = 77 t ( 16 ) 




\Jp 2 + mf 



are the dynamical equations for the 2-body system coupled to gravity, where 
% = 1,2. 

The defining equation for the Hamiltonian is, using eqs. flUD and (|i"5j), 



where B 1>2 = H - + mf >2 , J 2 = (y 1 H 2 + 8A/k 2 - 2epf - 8A/k 2 and 
r = Z\ — z 2 . For a given A > —(kH) 2 /8, equation (|T7|) describes the surface in 
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(r,p,H) space of all allowed phase-space trajectories. Since if is a constant 
of the motion, (a fact easily verified by differentiation of flT7| ) with respect 
to t) a given trajectory in the (r,p) plane is uniquely determined by setting 
H = H in (0). 



In the equal mass case the canonical equations of motion (|16|) can be 
solved exactly. The proper time of each particle is 



dr = dt N Q (z a) 



rn 



dt 



WkKqKxX 



m 



y/p 2 + m 2 J + KMi \/p 2 + m 2 

yielding upon insertion in (|T6|) the solution p = ^p(/(r) — l//(r)), where 



(18) 



f(r) 



" + —[J- V - — lm tan 

1 + 7™ 2 — — 



-~lm{T— tq) 



-7m(T-T ) 



H 
H 



to) 



7m < 



7n 



7m > 



(19) 



8A 



o- = (1 + a/7h)(vPq + m 2 -ep )-^ r , and 77 



2T7=' 



where 7z — It ^2^2? <-> — < • v " M V j " n "' ' '"" ' // 
with p Q the initial momentum at r = r . It is then straightforward to obtain 
an exact expression for r as a function of r, either directly from the second 
of eqs. (p~6|) or by inserting the expression for p(r) into flT7|) and solving for 
r. 

Analysis of these solutions shows that two types of motion are possible, 
depending on the amount of energy and the value of the cosmological con- 

2 4 

stant. For A < m A _ ^ = A c the two particles will always execute periodic 



2(_ff 2 -4m 2 ) 

motion, with period 



16 



T 



VK 2 m' 2 +8X 
16V H 2 -4m 2 
KHm 
16 

V-K 2 m 2 -8A 



tanh -1 ( ^ K ' 2m ' 2 + 8A ^ H2 ~^ m2 

\ kHtti 



tan 



— 1 / V — n 2 m 2 — 8AVH 2 -4m 2 



7m > 
7m = 
7m < 



(20) 



which can be obtained from the initial and final values of p when r = from 
(IIT]). If A < this is the only state of motion possible. However if A > 
there also exist a countably infinite set of unbound states of motion, in which 
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the particles begin with some (sufficiently large) separation, approach one 
another at some minimal value of |r|, and then move apart from one another 
toward infinite separation. 

For A > A c , only unbounded motion is possible. If the initial momentum 
is sufficiently small, the particles will cross one another before receding to- 
ward infinity. Otherwise they simply approach one another at some minimal 
value of |r| and then reverse direction toward infinity as previously men- 
tioned. One peculiar feature in this regime is that the two particles diverge 
to infinite separation at finite r! Specifically, if the above condition is sat- 
isfied, then the determining equation requires that r — > oo as p approaches 
the value ^ (V k 2 H 2 + 8A + v^A) . The corresponding value of r at which 
this will occur can be found from the 7 m > solution above. The explicit 
expression for this r is cumbersome and is omitted. 

In figure 1 we plot r vs. r for two bodies initially at r = for fixed 
A = —1.5 and H = 16 for several different values of m. As the motion 
becomes more relativistic (i.e. m gets smaller) we find that a second maxi- 
mum develops in the curve, which then vanishes for very small m. We find 
this behaviour for a broad range of values of A < and H. Figure 2 con- 
tains an analogous plot for A = 1.5, showing the transition from bounded 
to unbounded motion. As A — > A c (i.e. as m becomes small) the particles 
rapidly separate, remaining nearly stationary for an increasingly large period 
of proper time before coming together again. At A = A c this separation time 
becomes infinite, and for A > A c , the separation diverges at finite r. 

To summarize, we have obtained an exact solution to the 2-body problem 
in relativistic lineal gravity in which spacetime is not flat, in contrast to other 
lower-dimensional solutions to this problem in which spacetime is either flat 
outside matter || [7|, H or the matter sources are stationary || . The system 
is described by a conservative Hamiltonian in the canonical formalism. Many 
interesting features of this problem remain to be explored, including an ex- 
tension to the unequal-mass case, the transition from bounded to unbounded 
motion for A > 0, the appearance/disappearance of a second maximum for 
A < 0, and an investigation of the statistical properties of the model for large 
N. In this last case many-body forces play an important role - as shown in 
ref . Q , in the post-Newtonian approximation n-body forces appear at order 
K n_1 . We intend to make these the subjects of future investigation. 

This work was supported in part by the Natural Sciences and Engineering 
Research Council of Canada. 
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